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Abstract. 

We present a model for describing nuclear matter at finite density based on quarks interacting 
with chiral fields, a and vr. The chiral Lagrangian also includes a logarithmic potential, 
associated with the breaking of scale invariance. We provide results for the soliton in vacuum 
and at finite density, using the Wigner-Seitz approximation. We show that the model can reach 
higher densities respect to the Linear-a model, up to p « 3/9o for rria = 1200 MeV. 
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1. Introduction 

^ . The phase diagram of nuclear matter can be in principle more complicated than the one based 

0^ \ on MIT-bag-like models which predict a direct transition from hadronic matter to quark-gluon 

^^ ■ plasma (QGP). In the region of high densities and relatively low temperatures a variety of phases 

can exist in which chiral-symmetry breaking is realized in different ways. On the other hand, 
J ^ to study systems at finite density by using chiral lagrangians is not a trivial task: for instance 

(^ ' models based on the Linear o"- model fail to describe nuclear matter already at /> ~ po- In Ref. [1] 

the authors conclude that the failure of the a model is due to the restrictions on the scalar field 
dynamics imposed by the Mexican hat potential. In Ref. [2] they use a non-linear realization 
of chiral symmetry in which a scalar-isoscalar effective field is introduced, as a chiral singlet, to 
simulate intermediate range attraction. In this way the dynamics of the chiral singlet field is 
;J] ■ no more regulated by the Mexican hat potential and accurate results for finite nuclei can be found. 

Another possible solution to the problem of studying systems at finite density in chiral models 
is to use a linear realization but with a new potential, which includes terms not present in the 
Mexican hat potential. A possible guideline in building such a potential is scale invariance. 

In QCD scale invariance is spontaneously broken due to the presence of the parameter Kqcd 
coming from the renormalization process. Formally, the non conservation of the dilatation current 
is strictly connected to a not vanishing gluon condensate [3], [1], [5], [£]: 

^d^^QCo) = ^{F^^A^)F^"'{x))- (1) 

In the approach of Schechter, Migdal, and Shifman [7] a scalar field representing the gluon 
condensate is introduced and its dynamics is regulated by a potential chosen so that it reproduces 
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(at mean-field level) the divergence of the scale current that in QCD is given by Eq. ([T]). The 
potential of the dilaton field is therefore determined by the equation: 

e;^ = W{ct>) -<l^^ = 4.e,ac (^\ (2) 



where the parameter e^ac represents the vacuum energy. To take into account massless quarks a 
generalization was proposed in Ref . [8] , so that also chiral fields contribute to the trace anomaly. 
In this way the single scalar field of Eq. ([2]) is replaced by a set of scalar fields {o", tt, (/>}. 

It has already been shown that an hadronic model based on this dynamics provides a good 
description of nuclear physics at densities about pQ and descibes the gradual restoration of chiral 
symmetry at higher densities [9]. 

The new idea we develop in this work is to interpret the fermions as quarks, to build the 
hadrons as solitonic solutions of the fields equations and finally to explore the properties of the 
soliton at finite density using the Wigner-Seitz approximation. 

The structure of the paper is the following. In Sec. [2] we describe the model we are using, 
in Sec. [3] we discuss the results for the single soliton, first in vacuum and then at finite density. 
Finally, in Sec. H] we present our conclusions and future outlooks. 

2. The model 

The basic approach of a hybrid chiral model is to couple quarks to mesons in a chirally invariant 
way. We consider such a Lagrangian: 

C = i>{ird^, - gA^ + i^ ■ -ri^))^ + \{^^.<y^^a + d^n ■ d^^iz) - y(0o,^,vr). (3) 

Here if) is the quark field, a and vr are the chiral fields and (j) is the dilaton field which, in the 
present calculation, is kept frozen at its vacuum value (J)q. The potential is given by: 

y((/>0,o-,7r) : 
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where the logarithmic term generates from ([2]). The constants in the last line of Eq. (jj]) ensure 
that the vacuum energy is zero. 

The constants B and (J)q can be fixed by choosing value for the mass of the glueball and the 
vacuum energy e^^c) while 5 = 4/33 is provided by the QCD beta function and it corresponds to 
the relative weight of the fermionic and of the gluonic degrees of freedom. 

The first two terms of the potential are responsible for the breaking of scale invariance, while 
the term in the third line explicitly breaks the chiral invariance of the lagrangian. The Euler- 
Lagrangian equations that follow from Lagrangian ^ are: 



[i(^ - gT,{a + i-TT ■ tj^)]^' = 0, 
dV 
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The values of the parameters used for calculations in vacuum are listed in Table [TJ 



Table 1: Values of the parameters. 



Quantity 


Value 


|e^a,|V4 (MeV) 


236 


g 


5 


m^ (MeV) 


550 


B 


24.37 


e^^ (MeV) 


114 


6 


4/33 


(To (MeV) 


93 


(t)o (MeV) 


175.23 



2.1. The hedgehog ansatz 

We are working in the mean-field approximation, where mesons are described by time- 
independent, classical fields and where powers and products of these fields are replaced by powers 
and products of their expectation values. The quark spinor in the spin-isospin space is: 
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Once fixed the form of the quark wave function the self-consistents ansatze for the meson fields 
are: 

{&) = a{r), (wa) = ra7r(r) 

where cr{r), 7r(r) are radial functions of r. 



3. Results 

3.1. Single soliton in vacuum 

Solving the vacuum case requires the following boundary conditions for the fields: 



n'(0) = v{0) = 0, 
a'{0) = 7r(0) = 



(7) 
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Figure 1: Fields in vacuum for the solution g = 5, rricr = 550 MeV. 



while at infinity (in practice at a value r = R = A fm), the boundary conditions read: 

<T(i?)=ao,^(i?)=0, 



v{R) _ -ga{R)+e 



u{R) 



-ga{R)-e' 



(8) 



A test for convergence of the solution comes from another way of expressing the energy obtained 
by Rafelski |1U) by integrating out the fermionic fields: 
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Our solutions satisfy this consistency test up to a precision of the order of 10 ^. 



In Table [2] we present the static properties of the nucleon at Mean Field level and we compare 
them with experimental values and in Table [3] we show the decomposition of the soliton total 
energy in various contributions and the comparison with the Linear-o" model [11) . 



Table 2: Various nucleon properties at mean field level in the present work and comparison with experimental 
values. 



Quantity 


MFA 


Exp. 


M (MeV) 


1175.6 


1085 


{r!)i=o 


(0.73 /m)2 


(0.72 /m)2 


{rm)l=i 


(0.91 /m)2 


(0.80 /m)2 


9A 


1.27 


1.26 



Table 3: Contributions to the soliton total energy in the Logarithmic model and comparison with Linear-cr 
model. All quantities in MeV. 



Quantity 


Log. Model 


Linear-cr Model 


Quark eigenvalue 


83.1 


107.4 


Quark kinetic energy 


1138 


1056.9 


-Bct (mass+kin.) 


334.5 


320.3 


Etj (mass+kin.) 


486 


373.1 


Potential energy a — tt 


105.7 


120.7 


Eqa 


-101.4 


-62.3 


Eqw 


-787 


-673.2 


Total energy 


1175.6 


1136.2 



3.2. Single soliton at finite density 

In order to describe a soliton system at finite density we use the Wigner-Seitz approximation. 

This approach is very common for this kind of calculations and it has already been widely applied 

to both non- linear |12| [T3l [T3] and Linear-u models |15) . 

Specifically the Wigner-Seitz approximation consists of replacing the cubic lattice by a spherical 

symmetric one where each soliton sits on a spherical cell of radius R with specific boundary 

conditions imposed on fields at the surface of the sphere. The specific configuration of the meson 

fields, each one centered at a lattice point, generates a periodic potential in which the quarks 

move. 

In particular now the spinor of quark fields must satisfy the Bloch theorem: 
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(10) 



where k is the crystal momentum (for the ground state is equal to zero) and $fc(r) is a spinor 

that has the same periodicity of the lattice. 

We should explain in detail our choice of boundary conditions, because there are various sets 

of possible boundary conditions |15|, [H] . In particular we relate the choice of our boundary 

conditions to the parity r — )• —r operation; respect to this symmetry the lower component v{r) 

of quark spinor, the pion 7r(r) and the p{r) fields are odd, so all of these fields should vanish at 

R: 

v{R) = tt{R) = p{R) = 0. (11) 

For the remaining fields and the upper Dirac component we apply the usual "fiatness" condition 
on the boundary: 

u'{R) = a'{R) = J{R) = A'siR) = A't{R) = 0. (12) 

Basically the calculation is based on solving the set of coupled fields equations in a self- 
consistent way for each step in R; we start from the vacuum value, R = 4 fm, and we slowly 
decrease the cell radius down to the smallest radius for which self-consistent solutions can be 
obtained. 
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Figure 2: Upper and lower components of the Dirac spinor, a and pion fields for different values of the cell 
radius R. 

In Fig. [2] we plot the Dirac and the chiral fields for different values of i?; up to i? = 2 fm, the 
solutions do not change significantly, but as the cell radius shrinks to lower values, we see that 
all the fields are deeply modified by the finite density. In Fig. [3] we show the results, for the total 
energy of the soliton, in the present model and in the Linear-o" model. For a fixed value of the a 
mass, the logarithmic model do reach higher densities; as the a mass raises, the system can get 
to lower K because chiral fields are stuck on the chiral circle. 
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Figure 3: Total energy of the soliton as a function of cell radius R for the Linear-cr model [15] for different values 
of ttIct and present model. 



4. Conclusions 

We used a Lagrangian with quarks degrees of freedom based on chiral and scale invariance to 
study how the soliton behaves in vacuum and at finite density. At zero density the interplay 
between quarks and chiral fields lead to a soliton of mass M = 1175.6Mey, not too far from 
the experimental value given by the medium value 1085 MeV of M^ + M/\. We employed the 
Wigner-Seitz approximation to describe the dense system and we compared our results to the 
Linear-o" model |15| . In this approach we showed that the new potential, which includes the scale 
invariance, allows the system to reach higher densities at fixed m„. This improvement is totally 
due to the different dynamics between chiral fields, given by the logarithmic potential. 

The present work will be extended in several directions. First of all the finite density approach 
will be applied to a model including vector mesons, in order provide the necessary repulsion at 
large density. Next we will provide a precise and accurate calculation of the eigenvalues band in 
the soliton crystal in order to see how the system is affected by these effects |15) . Finally, the 
model can also be studied also at finite temperature, including the dynamics of the dilaton field. 
It is interesting to note that when this model has been investigated at finite temperature assuming 
the fermions to be hadrons [9], a phase diagram similar to the one proposed by McLerran and 
Pisarski |16| was obtained. In principle the approach based on the Wigner-Seitz scheme should be 
able to recover a scenario similar to the one discussed in [9] , but starting from more fundamental 
ingredients. 
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